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Abstract. Self-organised criticality (SOC) has been sug- 
gested as a potentially powerful unifying paradigm for 
interpreting the structure of, and signals from, accretion 
systems. After reviewing the most promising sites where 
SOC might be observable, we consider the theoretical ar- 
guments for supposing that SOC can occur in accretion 
discs. Perhaps the most rigorous evidence is provided by 
numerical modelling of energy dissipation due to magne- 
tohydrodynamic turbulence in accretion discs by G Geert- 
sema & A Achterberg (A&A 255, 427 (1992)); we inves- 
tigate how "sandpile"-type dynamics arise in this model. 
It is concluded that the potential sites for SOC in accre- 
tion systems are numerous and observationally accessible, 
and that theoretical support for the possible occurrence 
of SOC can be derived from first principles. 

Key words: accretion; accretion discs - stars: cataclysmic 
variables - galaxies: Seyfert - X-rays: stars 



1. Introduction 

The question whether accretion discs can be in a state of 
self-o rgani sed criticality (SOC) (Bak ct al. 1988 , Kadanoff 
et al. 1989 ) was raised explicitly by Mincshige ct al. ( 1994 ), 
and strong ly implicitly by Bak et al. ( 1988 ). Young & Scar- 
gle (1996) have raised the related question of transient 
chaos in accretion systems. A distinctive feature of SOC 
is flickering energy transport with no characteristic length- 
scale or time separation, displaying a 1// power spectrum; 
in this context, 1/ / is shorthand for inverse power law fre- 
quency dependence with unspecified index. In SOC sys- 



tems, for which mathematical sandpiles (Bak et al. 1988, 
Kadanoff et al. 1989) provide a paradigm, global transport 
occurs as a result of self-organised avalanches which are 
triggered locally when the accretion of sand leads to a crit- 
ical gradient being exceeded at a given point. This causes 
local redistribution of sand, which may lead to the critical 
gradient being exceeded at neighbouring points, result- 
ing in further redistribution and, cumulatively, to a global 



avalanche. Following the avalanche, the system returns to 
a subcritical configuration; accretion then continues until 
it again creates a local excess gradient, triggering a further 
avalanche. In the SOC-sandpile paradigm, it is helpful to 
note that the word "critical" is being used in two senses: 
the sandpile has a critical gradient, redistribution being 
triggered wherever this gradient is exceeded; and the self- 
organised global avalanches which emerge from the inte- 
grated effects of local redistribution may be scale-free, in 
which case they are linked to the generic field of critical 
phenomena. 

In mathematica l mod els and also some experi mental 
realisations (Nagel 



1992 



1995 



Frette et al. 



1996 



Feder 

Christensen et al. 1996 ), avalanche statistics display scale- 
free 1/f characteristics. SOC sandpile algorithms are ex- 
tremely simple and possess the attraction of any success- 
ful reduced system: it becomes unnecessary to attempt 
to model the detailed, and perhaps insuperably complex, 
microphysics of transport in the real system. Such an ap- 
proach has been applie d to the statistics of solar X-ray 
bursts (Lu 



1991 



for a recent treatment, see 



Hamilton 

MacKinnon & Macpherson |1997| and references therein 



) and to transport phenomena in magnetically confined 
plasmas (Newman et al. 1996| , Carreras et al. 1996, Dendy 



k Helander (19971 
lish whether 



1998) 



It is clearly of interest to estab- 
in certain circumstances, the detailed mod- 



elling of astrophysical accretion flows could be substituted 
by a simple SOC sandpile paradigm. This is a highly 
cross-disciplinary question. In the present paper, we aim 
to carry forward the debate in two ways: by identifying 
in greater detail the classes of accretion flow where SOC 
might play a role; and by examining theoretical arguments 
that we believe point clearly towards SOC in some accre- 
tion flows. 



2. Some sites of interest 
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It was noted by Bak et al. ( |1988| ) that flickering signals 
with 1/f power spectra have been observed for the X- 
ray variability of active galactic nuclei (AGNs), specifi- 
cally: 0.05-2keV X-ra ys from the Seyfert galaxy NGC4051 
(Lawrence et al 1987) and 2-7keV X-rays from the Seyfert 
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galaxy NGC5506 (McHardy & Czerny |1987| ), and for 10- 
140 keV X-rays from the massive compact binary Cyg X-l 
(Nolan et al. 1981). Similar 1// spectra for X-ray vari- 



ability from binary accreting systems were noted by Mi- 



neshige et al. (1994) (neutron stars, Makashima 1988) and 
by Geertsema & Achterberg (1992) (cataclysmic variables 
and dwarf novae, Wade & Ward 1985). While we believe 



this list can be extended, as we discuss below, it is also im- 
portant to be more specific about where, in these diverse 
systems, SOC might be occurring. 

Perhaps the simplest case is presented by a paradig- 
matic AGN. Let us take for this the standard picture of a 
massive black hole fed by a cascade of structures: the ac- 
cretion disc, the molecular torus and on a larger scale the 
galactic disc and its barred or spiral structure, which to- 
gether generate a clumpy and irregular transfer of gas and 
stars. This irregular mass transfer could be analogous to 
the random and discrete feeding of a sandpile with grains 
of sand near its apex. It is believed that, at least in thin 
discs, the residence time T res of matter in the disc is much 
longer than the free-fall time, Tff w (R 3 /GM) 1 / 2 at ra- 
dius R: 



T ff 



R 2 
otH 2 



where H <C R is the thickness of the disc, and a <C 1 is 
the well-known viscosity parameter. The condition r res 3> 
Tff, which is implicit in most accretion disc models, ap- 
pears to be necessary, at least in principle, for the validity 
of a sandpile-type approach to mass transfer within the 
disc. 

Thus if the accretion disc, like the sandpile, is in a 
state of SOC, there would be no clear link between the 
pattern of accretion from the torus and the pattern of 
avalanches leading to mass transfer across the disc, over 
its inner edge, and onto the black hole, giving rise to the 
X-ray signal. The latter would automatically display 1/f 
flicker. On larger space and timescales relating to our own 
Galaxy, even steady gas inflow from the Galactic bar could 
also result in avalanches to the inner regions, and these 
avalanches could give rise to episodes of AGN activity dur- 
ing an otherwise quiescent phase, with a duty cycle of a 
few percent, as inferred on statistical grounds for other 
AGNs (Mezger et al., |1996| ). The sandpile-SOC paradigm 
thus provides a candidate framework for summing up the 
consequences of the complex and intricate physics that re- 
lates the largescale dynamics of the Galactic disc to the 
activity of its central black hole. 

Cataclysmic variables, and in particular dwarf novae, 
present further opportunities for SOC. First, the flow of 
matter from the secondary across the inner Lagrangian 
point could itself be an avalanching process governed by 
SOC. In this case, if the radiation source was a hot spot 
where the accreting mass flow reached the outer edge of an 
accretion disc, it would flicker with 1/f statistics. Steady 
mass flow from the secondary would also be compatible 



with an SOC signal, however, provided the latter origi- 
nated from SOC mass transfers from the accretion disc 
(or accretion column in the strongly magnetised regime) 
to the white dwarf. Either regime appears possible in prin- 
ciple, both for dwarf novae and in the wider context of bi- 
nary accreting systems. In dwarf novae, it is pointed out 
by van Amerongen et al. ( |1990| ) and Lasota et al. (1995) 
that outbursts are due to suddenly increased accretion to- 
wards the white dwarf, but that it is unclear whether the 
instability resides in mass transfer from the secondary to 
the accretion disc, or across the accretion disc and onto 
the white dwarf. We note that, observationally, flicker- 
ing is more often strong during the quiescent phase of 
dwarf novae than during major eruptions (Wade & Ward 
1985 ). Questions relating to the locatio n of un stable flows 
in soft X-ray transients (Lasota et al. |1996| ) are similar 
to those in dwarf novae. For the wind-fed X-ray binary 
pulsar G X301-2 , a model has been proposed (Orlandini 
& Morrill 1992) involving "noisy" accretion of blobs of 



matter formed by magnetohydrodynamical (MHD) insta- 
bility at the magnetospheric radius, and not caused by 
inhomogeneities present in the stellar wind from the opti- 
cal companion. This approach is somewhat reminiscent of 



that of Baan (1977), where accreted matter accumulates 
at the magnetopause of a rotating neutron star until an in- 
terchange instability is triggered, after which the released 



matter generates an X-ray burst. Baan (1977) suggested 
that, since most of the time between bursts is a refilling 
time, an approximately linear relation should exist be- 
tween burst energy and the subsequent quiescent interval. 
However, in an SOC model where randomly arising local 
instability is sufficient to trigger a global avalanche, there 
would be no such correlation. This appears to be a key ob- 
servational discriminant for the possible presence of SOC 
in a given accreting system. We also n ote t hat, from the 
theoretical point of view, Frank et al. 1992 have pointed 



out that a local instability in a given annulus of the disc 
can only trigger large-scale instability accross the disc if 
parameters in neighbouring annuli are such that the ef- 
fect of this instability in these annuli can in turn trigger 
local instability there. This amounts to a prescription for 
a sandpile-type approach, and hence for the possibility of 
SOC. 

It seems clear from the foregoing that there is good 
observational and interpretative motivation for testing for 
SOC in a broad range of accreting astrophysical systems, 
encompassing flickering AGNs and certain distinct loca- 
tions within a variety of binary objects. Firm identifica- 
tion of SOC would yield information on the global conse- 
quences of the smaller-scale physics of the accretion pro- 
cess, while short-circuiting the need for detailed modelling. 
Let us now turn to theoretical arguments for expecting 
SOC. 
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3. Theoretical considerations 

The identification of mechanisms responsible for the dis- 
sipation of shear flow energy within an accretion disc is 
an active field of research, because of its importance in 
determining transport across the disc and onto the com- 
pact object. If, as has been suggested (Bak et al 



1988 



Mineshige et al. 1994), the combined global effects of lo- 
cal transport physics result in SOC, observational signa- 
tures of the type described above will emerge. It is widely 
accepted (see, for example, Longair (1994) and Narayan 
(1997)) that anomalous viscosity caused by MHD turbu- 
lence probably plays an important role in the flux of an- 
gular momentum within accretion discs. MHD turbulence 
can arise naturally in accretion discs, see for exam ple th e 
instability mech anism s proposed by Tagger et a l. (199lj ), 
Vishniac et al. (|l99C|) , and Balbus & Hawle y (|l99l|) . It 



has also been pointed out (Chen et al. 1995 ) that there 
exist ranges of accretion rate M and disc radius R for 
which no stable steady state solution of the basic equi- 
librium equa tions is possibl e. In t his case, it is suggested 
(Chen et al. 1995 , Narayan 1997 ) that the flow is forced 
into a time-dependent variable mode, satisfying the re- 
quired AI in the mean. This provides further motivation 
to test the applicability of sandpile-type models. In this 
section, we concentrate on the model for MHD turbulent 
energy dissipation in accretion discs presented by Geert- 
sema & Achterberg (1992), because it appears to give rise 
to SOC. While the link to cataclysmic variable and dwarf 
nova observations (but not, explicitly, SOC) was made by 
Geertsema & Achterberg (1992) themselves, we believe 
that this work has much wider implications: for the gen- 
eral question of SOC in accretion flows; for the connec- 
tion between SOC and turbulence in general; for the role 
of SOC in plasma physics, since it represents the first in- 
stance where SOC has been observed in a mathematical 
model derivable from the fundamental equations of MHD; 
and for the role of SOC in tcrrestial experimental systems 
- real sandpiles and ricepiles, as distinct from mathemat- 
ical idealisations thereof - where uncertainties about its 



scope remain (Na gel [1992] , Feder |l995i Frette et al. |1996 , 
Christensen et al. 1996 ). 

Before considering the model of Geertsema & Achter- 
berg ( 1992 ) in greater detail, let us turn t o its results. Fig- 
ure 12 of Geertsema & Achterberg (1992) shows the calcu- 
lated times series of energy dissipation events within the 
disc. We note that this is qualitatively very similar to the 
observed time series of energy dissipation measured in an 
experimental ricepile displaying SOC, Fig 2c of Frette et 



al. (1996), and in a related mathematical model of Dendy 
& Helander ( |1997| , fig. 3, and |1998| , fig.5) More quantita- 
tively, Fig 13 of Geertsema & Achterberg (1992) shows the 



power spectrum of energy dissipated by MHD disc turbu- 
lence, which displays the 1// dependence characteri stic o f 
SOC; compare, for example, Fig 3 of F rette et al. ( 1996 ) 
and Fig 3 of Christensen et al. ( 1996 ) which show mea- 



sured spectra of energy dissipation and particle transit 
times, respectively, in SOC ricepiles. 

Given the clear indications of SOC emerging from 
the MHD turbulence model of Geertsema & Achterberg 



(1992), it will be of interest to establish how it has arisen. 
A full explanation must await diagnostic analysis of the 
code runs generated by this model. Pending this, we con- 
clude the present section by seeking to identify some of the 
relevant salient features. In outline, the model of Geert- 
sema & Achterberg ( 1992| ) is constructed as follows. 

The accretion disc is regarded as a differentially rotat- 
ing turbulent MHD fluid, and is modelled by a reduced 
system of equations reflecting the most important fea- 
tures of three-dimensional MHD turbulence. The disc is 
assumed to be thin in comparison with its diameter, and 
the flow is taken to be subsonic and hence incompress- 
ible. The possible existence of a large-scale magnetic field 
is neglected, and the only turbulent structures considered 
have a length scale shorter than the height of the disc, 
making the turbulence essentially three-dimensional. In a 
coordinate system rotating with the disc angular velocity 
= b z {Ik{R) at radius r = R, the simplified MHD equa- 
tions for the flow velocity u and the normalised magnetic 
field vector b = B/(47rp) 1 / 2 are 



-Vn-2J7 x u- 
+ vV 2 u 



K 



(It 



R 



b • Vb - u • Vu, 



K 



dt 



b X Gy 



rjW 2 b + b • Vu u ■ Vb, 



where x = r — R and y — Rcj) are the local radial and 
azimuthal coordinates respectively, rj is the resistivity, and 
v is the viscosity coefficient, n = P/p + 6 2 /2 is the total 
scalar pressure, while the b ■ Vb term contains the off- 
diagonal terms of the MHD stress tensor. 

As these equations are still too complex for a numer- 
ical treatment over a sufficiently large dynamic range in 
scales, they were simplified further following a suggestion 
by Desnyanski & Novikov (1974) in the context of hydro- 
dynamics, and later applied to MHD by Gloaguen et al. 



(1985). In this approximation, the space of wave vectors 
k is discrctized into a finite set of k n , and the non-linear 
interaction between different components of the Fourier 
transforms of u and b is described by the set of equations 



du n 
~~dT 



vk n u n 



= ak n (i4_! - 2u n u n+1 ) 

+ fik n (u n -iUn - 2w 2 l+1 ) - (u <-> b), 



— ^ + r\k 2 n b n = ak n (m„_i&„_i - 2u n b n+1 ) 
dt 

+ (3k n (u n -ib n - 2u n+1 b n+ i) - (b<-+u). 

In the three-dimensional generalisation of this approxi- 
mation scheme, the discretisation is made by dividing the 
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k-space into into spherical shells, thus discarding the infor- 
mation regarding the direction of k. This allows a greatly 
simplified system of nonlinear equations to be written 
down, analogous to that of Gloaguen et al. (1985), which 



are supplemented with additional terms to account for the 
effects of differential rotation. 

These equations, which are taken to model the turbu- 
lent cascade of MHD, were solved numerically by Geert- 
sema & Achterberg (1992). They found that the turbulent 
shear stress can be very large, and has large, chaotic fluc- 
tuations on time scales of a few rotation periods. Perhaps 
the most striking feature of the simulations is, however, 
that the dissipation of energy at the smallest scales of the 
turbulent cascade is very intermittent. Energy is released 
in avalanches with a wide range of sizes, and the power 
spectrum of the dissipation rate obeys a 1// power law 
over nearly two orders of magnitude in intensity, see their 
Fig 13. This behaviour is, as already stated, similar to that 
of simple mathematical sandpile models. While a power 
law is to be expected from any scale-free model, in partic- 
ular a Kolmogorov-type one in the inertial range, we note 
that the similarities between this system and mathemat- 
ical sandpiles apparently extend further. Both are funda- 
mentally governed by nearest-neighbour interactions be- 
tween a discrete number of nodes. In the MHD accretion 
model, these reside in k-space, so the interaction is be- 
tween wave modes with similar wavelengths rather than 
between adjacent regions in real space. The dissipation is 
provided by viscosity and resistivity at large k, resembling 
the removal of material from the edge of a sandpile. The 
MHD accretion model is, of course, much more complex 
than the simple sandpile algorithms considered so far in 
the literature. Two fields, u and b, each with three compo- 
nents, are involved rather than the single height parameter 
of conventional sandpile algorithms, and the time evolu- 
tion is governed by differential equations rather than dif- 
ference equations. Nevertheless, both the MHD accretion 
disc and simple sandpile models appear to exhibit similar 
self-organised, critical behaviour, supporting the claim of- 
ten made in the sandpile literature that SOC is universal 
phenomenon shared by large classes of cellular automata. 

4. Discussion 

We have shown in this paper that the numerical simu- 
lations by Geertsema & Achterberg (1992) of viscous re- 



sistive MHD turbulence in an accretion disc give rise to 
behaviour characteristic of self-organised criticality in a 
sandpile. This similarity may help to explain certain ob- 
served properties of a range of accreting astrophysical sys- 
tems, which we have reviewed. Furthermore, the result is 
of intrinsic scientific interest as an ab initio demonstra- 
tion of the emergence of SOC from a system of MHD- 
based equations. Further analysis of the numerical re- 
sults, to establish how the SOC-sandpile phenomenology 
arises, would be of great interest both for accretion disc 



astrophysics and for fundamental plasma physics. At the 
present stage of development of realistic MHD simulations 
(for example of the link between smallscale dissipation and 
the time behaviour of the flux of matter onto the central 
object), direct numerical proof of the occurence of SOC- 
sandpile phenomenology is still a remote objective. How- 
ever, the capacity of the SOC-sandpile paradigm to cir- 
cumvent complex "full mathematical descriptions" is the 
basis for its present attraction in many fields of physics. 
If future work confirms the importance of SOC-sandpile 
phenomenology in accreting systems that is suggested in 
the present paper, it would permit a dissociation (at least 
at zeroth order) of the detailed physics of turbulence in 
the disc from the global modelling of this class of astro- 
physical object. 
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